Prediction of superconducting properties of CaB2 using anisotropic Eliashberg theory 
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Superconducting properties of liypothetical simple hexagonal CaB2 are studied using the fully 
anisotropic Eliashberg formalism based on electronic and phononic structures and electron-phonon 
interactions which are obtained from ab initio pseudopotential density functional calculations. The 
superconducting transition temperature Tc, the superconducting energy gap A(k) on the Fermi 
surface, and the specific heat are obtained and compared with corresponding properties of MgB2. 
Our results suggest that CaB2 will have a higher Tc and a stronger two-gap nature, with a larger 
A(k) in the a bands but a smaller A(k) in the tt bands than MgB2. 
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The discovery of superconductivity in MgB2" with its 
remarkable superconducting transition temperature 
as high as 39 K has triggered many investigations of 
related materials including various binaries: BeB2r^ 
CaB2,^^^>^i NaB2,-^ SrB2,-i>i CuB2,-^ ZrB2,^i^ TaB2,-^^ 
OsB2,^ ScB2,-^'i YB2^ MgBe2,-^ MgB4,-^ and LiBi^; 
and ternaries: CaBeSi^ UBC^^^^^^^^ NaBC^ 
MgB2_,C„^^ Mgi_,Ca,B2f^'^'il^ Mgi_,Li,B2f^'^ 
Mgi_,Na,B2^'^ MgBe,B2_,^^ CuB2_,C,^ MgB2C2^ 
and Be2Ba;Ci_a;ji^ For these compounds, their 
structural, electronic, and superconducting prop- 
erties are studied experimetally^ i^^'^^i^^'^^i"'^^ and 
theoretically^^^^^^i^^^^^^^ia The structural 
stabilities of the existing and hypothetical compounds 
are studied extensively as weWJ^ Even with these 
efforts, however, the highest Tc observed experimentally 
in these compounds is still that of MgB2. There have 
been theoretical predictions of higher Tc's than 39 K, for 
example, - 50 K in CuBa-xCj and Tc - 100 
and 65 Kii in hole-doped LiBC. 

One of the interesting candidates for a higher Tc is 
the hypothetical simple hexagonal CaB2, with the same 
crystal symmetry as MgB2, with all Mg atoms being re- 
placed by Ca atoms. First-principles calculations predict 
that the simple hexagonal CaB2 should have a greater 
unit cell (both a and c lattice constants are longer) and a 
much larger value for the density of states (DOS) at the 
Fermi energy (Ep) than MgB2.-'^'^'^'^ Both features are 
favorable for a higher Tc, since MgB2 is found to have a 
positive dependence of Tc on the unit-cell volume^! and 
more generally, a larger DOS at Ep can result in a greater 
electron-phonon coupling constant. Although CaB2 is 
an attractive candidate for superconductivity, it has not 
been synthesized as yet, nor has it been extensively stud- 
ied theoretically for superconductivity. Since the Fermi 
surface of CaB2 is predicted to consist of multiple sheets 
similar to those in MgB2^i^'^i^ it will be necessary to use 
the anisotropic Eliashberg theory to deal with the effects 
of possible variations of the superconducting energy gap 
on the Fermi surface as in the case of MgB2;22,:i? 

In this work, we perform density-functional ab-initio 




FIG. 1: Atomic structure of hypothetical simple hexagonal 
CaB2. The simple hexagonal unit cell contains one Ca and 
two B atoms. 



pseudopotential calculations to determine the lattice con- 
stants of the hypothetical simple hexagonal phase of 
CaB2, and we calculate its electronic structure, phonon 
spectrum, and electron-phonon properties. From these 
material properties, we construct the fully anisotropic 
Eliashberg equations to obtain the superconducting prop- 
erties such as Tc, the momentum-dependent supercon- 
ducting energy gap A(k, T), and the specific heat as func- 
tions of temperature T. These results predict that CaB2 
has a higher Tc than MgB2 and a larger difference be- 
tween the superconducting energy gap in the a and tt 
bands when compared with MgB2. 

We consider the hypothetical simple hexagonal CaB2 
shown in Fig. 1. The atomic structure has exactly 
the same symmetry as MgB2, having one Ca and two 
B atoms in a unit cell. We determine the lattice 
constants a and c by minimizing the total energy of 
the system obtained from first-principles pseudopoten- 
tial density-functional calculationsj^i^ In the calcula- 
tions, plane waves of energy up to 60 Ry are used to 
expand the electronic wavefunctions, norm-conserving 
pseudopotentials^^ are used to describe the electron- 
ion interaction, and the local density approximation 
is employed to deal with the electron-electron interac- 
tions effectively. In more detail, we use the partial- 



FIG. 2; (Color online.) Electronic structure of CaB2: (a) 
electronic band structure and (b) the Fermi surface. In (b), 
the cylindrical (red and orange) sheets along the F— A— F line 
are of hole type, and the connected (blue) sheet along H-L 
lines is of electron type. 



core correction^ for the Ca pseudopotential, Ceperley- 
Alder exchange-correlation energyS^ of the Perdew- 
Zunger parameterization,^^ a 12 x 12 x 12 /c-point grid 
to evaluate the self-consistent electron density, and a 
18 X 18 X 12 fc-point grid for Fermi-surface properties. 

For simple hexgonal CaB2, the calculated lattice con- 
stants at ambient pressure are a = 3.185 A and c = 4.060 
A. Thus the c/a ratio is 1.27. Compared with theoret- 
ical values for MgB2 (which are a = 3.07 A and c = 
3.57 A, with c/a — 1.16), a in CaB2 is 4% greater and 
c is 14% greater. To check the ratio c/a at high pres- 
sure, we obtain the lattice constants at 30 GPa, which 
are a — 3.014 A and c = 3.842 A. Thus, the c/ a ra- 
tio is still 1.27 at this pressure. The obtained c/a — 
1.27 is larger than 1.165 known as the maximal c/a ra- 
tio of existing diboridesj^S and this is in agreement with 
previous calculationsi^'^i^'i Since the hexagonal CaB2 
has not been synthesized successfully, its formation en- 
ergy and stability have been studied theoretically^iii^S 
Hexagonal CaB2 is predicted to be stable with respect to 
hep Ca and rhombohedral B (a-Bi2)r but unstable to a 
phase separation into fee Ca and CaBg^ and less stable 
than (5-CaB2 phasci^ As reported theoretically, hexago- 
nal CaB2 may not be the lowest phase energetically, but 
our frozen-phonon calculations at high symmetry points, 
which will be presented below, do not show any struc- 
tural instability from the simple hexagonal phase. 

With the optimized lattice constants, we perform first- 
principles electronic structure calculations. Figure 2 
shows the obtained band structure along the high sym- 
metry lines and the Fermi surface. Overall, the band 
structure of CaB2 is similar to MgB2. In Fig. 2(a), the 
a and the tt bands are from the boron sp^ and pz or- 
bitals, respectively. The full bandwidth of the a bands 
is reduced substantially, compared with MgB2, because 
of the larger B-B bond length. As in the case of MgB2, 
the Fermi energy in CaB2 is lower than the top of the a 
bands, producing two cylindrical hole-type sheets along 
the F-A line [red and orange regions in Fig. 2(b)]. While 
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FIG. 3: (Color online.) Phonon dispersion relations and 
electron-phonon interactions in CaB2. (a) Dashed (blue) lines 
are the phonon dispersions within the harmonic approxima- 
tion and solid (red) lines are those considering anharmonicity 
in the frozen-phonon calculations. As in the case of MgB2, the 
doubly degenerate -E29 mode along the F — ^ line is highly 
affected by anharmonicity. In the harmonic case, the fre- 
quency of the -E29 mode is even lower than the A2U mode, 
but the former becomes higher than the latter in the anhar- 
monic case, (b) Phonon density of states F{u}) in meV~^ and 
the Eliashberg function a^F(uj) including anharmonic effects 
on the phonon frequencies. Solids (red) lines are the phonon 
density of states and the Eliashberg function in CaB2 . Dashed 
(blue) lines are those in MgB2 drawn for comparison. 



the hole-type sheets are close to those in MgB2 qualita- 
tively, the electron-type sheet from the tt bands in CaB2 
is quite different from the one in MgB2. In MgB2, the 
TT bands are above Ep at the AI and K points; however, 
in CaB2, they are below Ep. Thus, CaB2 has only one 
electron-type Fermi surface from the tt bands [the blue 
sheet in Fig. 2(b)], missing the Fermi surface along the 
K-L line that is present in MgB2. Hence, as shown in 
Fig. 2(b), the Fermi surface in CaB2 consists of three 
sheets: two hole-type cylindral ones along F-A lines and 
one electron- type along H-L lines. The density of states 
at the Fermi level in CaB2 is 0.96 states/eV per formula 
unit, which is much larger than 0.69 states/eV for MgB2. 
The contributions of the three sheets of the Fermi surface 
to the DOS at Ep are 12%, 32%, and 56% from the inner 
and outer hole- type sheets and the electron- type sheet, 
respectively. Our results for the band structures are in 
good agreement with the results using the full-potential 
linear muffin-tin orbital (FP-LMTO) method.^^i^ 

To consider electron-phonon interactions, we calcu- 
late the phonon structures in CaB2 based on the frozen- 
phonon method. In the first step, phonon frequencies and 
polarizations are obtained at high symmetry points of the 
Brillouin zone by frozen-phonon methods, and then in- 
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terpolated to the full Brillouin zone by interpolating the 
dynamical matrix. In the frozen-phonon calculations, 
we extract quadratic and higher-order dependences of 
the total energy on atomic displacements, and calculate 
harmonic and anharmonic phonon frequencies with and 
without using the harmonic approximation. 

Figure 3(a) shows the calculated phonon dispersions 
along the high symmetry lines, with and without anhar- 
monicity. The phonon frequencies are calculated by con- 
structing the dynamical matrices at the high-symmetry 
points and interpolating them along the high-symmetry 
lines. In our frozen-phonon calculations, the increase of 
the total energy with boron displacements are highly 
anharmonic for E2g modes (the B-B bond stretching 
modes) at the F and A points. This is consistent with 
the reported large anharmonicity at F in CaB2.-3i We find 
that the E2g frequency is 45.2 mcV in the harmonic case 
and 55.8 meV in the anharmonic case, both of which are 
substantially lower than those in MgB2. All the phonon 
frequencies at the high symmetry points are positive, and 
therefore do not indicate structural instabilities. 

The upper panel in Fig. 3(b) shows the phonon density 
of states obtained by interpolating the dynamical ma- 
trix throughout the Brillouin zone with anharmonicity 
included. Compared with MgB2, the phonons in CaB2 
have lower frequencies. Phonon modes associated with 
Ca atoms should have lower frequencies because of heav- 
ier atomic mass of Ca relative to Mg. In addition, phonon 
modes associated with B atoms also have lower frequen- 
cies because of the elongation of B-B bonds. 

We calculate the momentum-dependent Eliashberg 
function Q!^i^(k, k', w) from the difference in the self- 
consistent potential with and without frozen phonons, 
using the anharmonic phonon frequencies, and then take 
the average on the Fermi surface to obtain the isotropic 
Eliashberg function a^Fiuj), as shown in the lower panel 
in Fig. 3(b). From a^F(a;), we can evaluate the average 
electron-phonon coupling constant X ~ 2 J duja'^F{uj)/uj 
and the logarithmically averaged phonon frequency = 
exp[(2/A) / duja^F{uj)\nuj/Lj]. This gives A = 0.69 and 
ioin = 50 meV, respectively. The value of A is 13% larger 
and LUin is 34% smaller than those for MgB2. 

We obtain the superconducting properties in CaB2 us- 
ing the fully anisotropic Eliashberg formalism. ^■^''^■^ The 
anisotropic Eliashberg equations at imaginary frequen- 
cies are given as 



Z{k,iuJn) = 1 + /„s„^Wk'A(k, k', n-n') 



v/c^^,+A(k',zc^„0' 
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FIG. 4: (Color online.) Superconducting properties in CaB2. 
(a) Superconducting energy gap A(k) on the Fermi surface 
at 4 K plotted using a color scale. The lower panel shows the 
distribution p( A) of the values of the energy gap as well as the 
color scale for the size of the energy gap. (b) Temperature 
dependence of the superconducting energy gaps [solid (red) 
lines]. Dashed (blue) lines are the superconducting energy 
gaps in MgB2 drawn for comparision. Vertical lines repre- 
sent the distributions of A(k) in CaB2 at 4, 20, and 40 K. 
Curves are fit to the separate averages, Aa and A^, where 
the average A-^ is greater than the peak energy of p(A) in (a) 
because of the high-energy tail of p(A) [green part in (a)], (c) 
The specific heat Ce over T of CaB2 [solid (red) line]. The 
dashed (blue) line is Ce/T of MgB2 drawn for comparision. 
Anharmonic effects on the phonon frequencies are considered 
in all calculations of the superconducting properties. 



where w„ = (2n + 1)ttT is the fermionic Matsubara 
frequency at temperature T, Z(]i,iLUn) and A(k, iw„) 
are the momentum-dependent renormalization function 
and the gap function, respectively, and A(k, k',n) repre- 
sents the momentum-dependent electron-phonon inter- 
action. Definitions of symbols and details of the nu- 
merical method are described in Refs. [H, [H, and [3^ . 
Here, we assume that the Coulomb pseudopotential fi* 
is isotropic, since its momentum dependence is not as 
strong as that of the electron-phonon interactionj^^i^i^l 
and assume ^*{ijJc) = 0.12 for the cut-off frequency ujc — 
0.5 eV, as assumed for MgB2<^'^ 

We obtain Tc by finding the highest T at which the 
anisotropic Eliashberg equations [Eqs. (1) and (2)] have 
a non-trivial solution for A(k). With the assumed value 
of /i*(wc) = 0.12, we obtain Tc = 48 K, which is about 
10 K higher than that of MgB2. To check the sensitivity 
of Tc to we also consider ij,*{uic) — 0.10 and 0.14, 
obtaining Tc = 50 and 46 K, respectively. Thus, we can 
expect a higher Tc in CaB2 than in MgB2 even in the 
case of relative large fi*. 

With Eqs. (1) and (2), and their analytic continuations 
to the real-frequency axis, we calculate the superconduct- 
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ing properties below Tc- Figure 4(a) shows the calculated 
superconducting energy gap on the Fermi surface at T = 
4 K, where the average values of A(k) on the three sheets 
of the Fermi surface are well separated from one another 
and they are 9.1, 8.2, and 1.5 meV, respectively, with 
the average on the two hole-type sheets being 8.4 meV. 
As shown in Fig. 4(b), when compared with MgB2, the 
average value of A(k) on the a bands (the hole- type 
sheets) is substantially larger but A^ for the tt bands is 
slightly smaller at low T, indicating an enhanced two-gap 
nature for CaB2. Figure 4(c) shows the calculated elec- 
tronic specific heat as a function of temperature. Because 
of the higher DOS at Ep and the larger A, the normal- 
state specific heat Cn above Tc is expected to be larger 
in CaBa [Cn = InT with 7„ = 3.82 mJmor^K-^) than 
in MgB2 (7,1 = 2.62 mJmol~^K^^). In addition, since 
Att in CaBa is smaller in magnitude than in MgBa, the 
specific heat at T < 10 K has a bigger hump in CaBa. 

To summarize, we obtained the lattice constants, the 
electronic and phononic structures, and electron-phonon 
interactions in hypothetical simple hexagonal CaBa using 



first-principles calculations. We then calculated the su- 
perconducting transition temperature Tc, the supercon- 
ducting energy gap A(k, T), and the specific heat using 
the fully anisotropic Eliashberg formalism. The obtained 
Fermi surface in CaBa consists of three sheets rather than 
four as in MgBa, and the phonon frequencies are lower 
than in MgBa. A substantially higher Tc is predicted for 
hexagonal CaBa, and furthermore the multiple supercon- 
ducting energy-gap feature is expected to be enhanced in 
CaB2, with larger A^ and smaller A^ than in MgBa. 
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